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ABSTRACT 

We investigate the star formation process in primordial environment in the presence of radiative 
feedback by other population III stars formed earlier. In this paper, we focus our attention on the 
effects by photodissociative radiation toward the full understanding of the radiative feedback effects. 
We perform three dimensional radiation hydrodynamics simulations on this issue as well as analytic 
estimates, paying special attention on the self-shielding effect and dynamics of the star-forming cloud. 
As a result, we find that the ignition timing of the source star is crucial. If the ignition is later than the 
epoch when the central density of the collapsing cloud exceeds ~ 10 3 — 10 4 cm -3 , the collapse cannot 
be reverted, even if the source star is located at < lOOpc. The uncertainty of the critical density comes 
from the variety of initial conditions of the collapsing cloud. We also find the analytic criteria for a 
cloud to collapse with given central density, temperature and the Lyman- Werner (LW) band flux which 
irradiates the cloud. Although we focus on the radiation from neighbor stars, this result can also be 
applied to the effects of diffuse LW radiation field, that is expected to be built up prior to the reionization 
of the universe. We find that self-gravitating clouds can easily self-shield from diffuse LW radiation and 
continue their collapse for densities larger than ~ 10 3 cm -3 . 

Subject headings: theory:early universe — galaxies: formation — radiative transfer — molecular 
processes — hydrodynamics 



1. INTRODUCTION 

One of the important objectives of cosmology today is to 
understand the way how the first generation stars or galax- 
ies are formed and how they affected the later structure 
formation, and how they reionize the surrounding media. 
In particular, the self-regulation of star formation in the 
first low mass halos is quite important, since it controls 
the Population III (POPIII) star formation activity in the 
early universe, which is the key for the reionization and 
the metal pollution of the universe. These problems have 
been studied intensively in the last decade, hence we al- 
ready have some knowledge on these issues (see review by 
Bromm & Larson 2004). However, studies which prop- 
erly address the radiation transfer effects are still at the 
beginning, in spite of their great importance. 

Radiation from the first generation stars play quite im- 
portant roles on the formation of stars/galaxies through 
two main physical processes. Since first sta r s are 
expected to be very massi ve jOmukai & Ni shi I [T998I: 
Abel, Bryan, Normanl l2000l: iBromm, Coppi fc Larson I 
20021 : iNakamura fc Umemural Il999l . 1200 if), the amount 
of ultraviolet photons are quite large. Once these kinds 
of stars are formed, 1) emitted Lyman- Werner (LW) band 
photons photodissociate the hydr ogen molecules in the 
neigh b orhood by Solomon proc es s (iHiaman, Rees fc Loeb I 

Brand I 120011 : 



feedbacks are basically expected to have negative effects on 
the structure formation (e.g. Susa & Umemura 2004ab). 
In contrast, enhanced fraction of electro ns by photoioniza- 
tion can lead to efficient H 2 forma tion (Shapiro & Kang 
119871: iKang fc Shapiro! Il992l : ISusa et all Il998l : 
lOh fc Haiman I |2002[ ) , and it could accelera t e the 
structure formation (Ricotti, Gnedin, & Shull 2001; 



19971: lOmukai fc Nishi I 



19991: iGlover 
l200lf ). Afterwards 



Machacek, Bryan fc Abel I2001D . Afterwards or simul- 
taneously 2) ionizing photons propagate into surrounding 
media, followed by the photoionization and heating of the 
gas. Since hydrogen molecules are the only coolant in low 
mass first generation objects, and the photoheated gas is 
too hot to be kept in such small dark halos, these radiative 



lO'shea, Abel Whalen fc Norman ll2005l:lNagakura fc Omukail 
| 2Q05t ISusa fc Umemural 120061 : lAhn fc Shapiro I 120061 : ~ 
lYoshida et al. Il2006f ). " 

Since these physical processes related to this issue are so 
complicated, we try to understand one by one. In this pa- 
per, we focus our attention on the photodissociation feed- 
back effect including the effect of self-shielding. We do 
not take into account the effects by phot oionization, which 
i s disc ussed with limited parameters in lSusa fc Umemural 
(|2006l ), and will be presented in detail in the forthcoming 
paper. 

Concerning the photodissociation feedback effects, there 
are two contexts. First one is the local feedback by 
the nearby POPIII stars formed in the same halo or in 
the neighbor halos. Another one is the feedback by dif- 
fuse LW background radiation, that is built up by the 
formed POPIII stars in many dark halos spread in the 
unive rse ([Machacek, Bryan fc Abel 1 l200ll : lYoshida et al. I 
120031 ). In reality, the former would particularly be im- 
portant at the early phase of the POPIII star formation 
episode, ahead of the latter effect. Local feedback might 
also play important roles in the star formation process in 
relatively larger halos whose virial temperature is higher 
than 10 4 K. The physical difference of the two cases are 
the different duration and intensity of radiation field. In 
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the latter case, intensity is not so high, whereas the ra- 
diation field irradiate the halos continuously since it is 
a background. In contrast, in the former case, the ra- 
diation source basically turn-off after the life-time of the 
source star, although the intensity of the radiation could 
be much higher than the diffuse background. In this paper, 
w e mainly investiga t e the local effects. 

lOmukai fc Nishi I (|l999h investigated on this issue in an 
uniform virialized halo. They argue that only a single mas- 
sive POPIII star can destroy H2 molecules within ~ lkpc, 
which is comparable to the virial radius of 1O 8 M ha- 
los at z > 10. This means that star formation in such 
low mass halo s is highly re gulated by themselves. On the 
other hand, iGlover fc Brand I (|200l[ ) considered more re- 
alistic clumpy halos. They simply assessed the photodis- 
sociation time scale and free-fall time of the dense clouds. 
Comparing these two time scales, they found that if the gas 
density is sufficiently high, dissociation time scale could be 
longer than the free-fall time, which allow the cloud to col- 
lapse. Although those arguments are quite clear, they are 
based upon very simplified configuration. Therefore, more 
sophisticated treatments are necessary to investigate this 
feedback process. In order to improve the understanding 
of this issue, we carry out three dimensional radiation hy- 
drodynamics simulations as well as more detailed analytic 
arguments. 

This paper is organized as follows: in the next section, 
we briefly summarize the numerical scheme employed. In 
setup of our numerical simulations are shown. We 
present the numerical and analytical results in g] and ^5] 
§6 and §7 are devoted to discussion and summary. 

2. NUMERICAL METHOD 

We perform numerical si mulations by Radiation- SPH 
code developed by ourselves (|Susa Il2006f ). In this section, 
we try to briefly summarize the scheme. The code is de- 
signed to investigate the formation and evolution of the 
first generation objects at z > 10, where the radiative 
feedback from various sources play important roles. The 
code can compute the fraction of chemical species e, H + , 
H, H - , H2, and by fully implicit time integration. It 
also can deal with multiple sources of ionizing radiation, 
as well as the radiation at Lyman- Werner band. 

Hydrodynamics is calculated by Smoothed Particle Hy- 
dro dynamics ( SPH) method. We use the version of SPH 
by lUmemural (1 19931) with the modification according to 
ISteinmetz fc Mullerl (1 19931). We adopt the particle resiz- 
ing formalism bv lThacker et al~l (|2000D . 

The non-equilibrium chemistry and radiative cooling for 
primordial gas are calcu lated by the code developed by 
ISusa fc Kitavamal (|2QQQh , where H2 cooling and reaction 
rates are mostly taken from iGalli fc Palla I (|1998| ) . As for 
the photoionization process, we e mploy so-called on-the- 
spot approximation (|Spitzerlll978D . but we do not get into 
the details since we do not show the results with ionizing 
radiation in this paper. 

The optical depth is integrated utilizing the neighbour 
lists of SPH partic l es. It i s similar to the code described in 
ISusa fc Umen^rTal ([2004a), but now we can deal with mul- 
tiple point sources. In our new scheme, we do not create 
so many grid poin t s on th e light ray as the previous code 
ISusa fc Umemural (|2004af ) does. Instead, we just create 



one grid point per one SPH particle in its neighbor. We 
find the 'upstream' particle for each SPH particle on its 
line of sight to the source. Then the optical depth from 
the source to the SPH particle is obtained by summing 
up the optical depth at the 'upstream' particle and the 
differential optical depth between the two particles. 

The opacity of LW radiation is c alculated by utiliz- 
ing t he original self-shielding function ([Draine fc Bertoldil 
Il996l ). The column density of H2 is evaluated by the 
method described above. We assume same absorption for 
the photons within LW band, i.e. between 11.26eV and 
13.6eV. 

The code is already parallelized with MPI li- 
brary. The computational domain is divid ed by so 
called Orthogonal Recursive Bisection method ( Dubinski~l 
Il996f ). The parallelization method for radiation trans- 
fer part is s imilar to Multiple Wave Fro nt method 

developed by Nakam oto, Umemura, fc Susa I (|200ll ) and 
iHieneman et al. I (|2005l ), but it is changed to fit the SPH 
scheme. The code also is able to hand le self-gravity 
with Barnes-Hut tree (|Barnes fc Hutl [l986), which is also 
parallelize d. The c ode is tested for various standard 
problems (Susal [20061 ) . We also take part in the code com- 
parison project w ith other radiation hydrodynamics codes 
(|lliev et aT1l2006[ ) in which we find reasonable agreements 
with each other. 



3. SETUP OF NUMERICAL SIMULATIONS 

We simulate the collapsing primordial gas cloud with 
M = 8.3 x 1O 4 M , which is initially uniform and spher- 
ical. The initial number density is 14cm -3 , whereas the 
corresponding radius is ~ 40pc. We employ four initial 
temperatures T ini = 50K, 100K , 200K and 350K which 
represent the coldness of the initial collapses. The uni- 
form cloud is embedded in an extended uniform low den- 
sity cloud with 0.1cm -3 and the same temperature as the 
dense cloud. Remark that such a thin surrounding mate- 
rial does not affect the fate of the dense cloud. The mass 
of an SPH particle is O.32M . We employ 524288 particles 
in total. Approximetely half of the particles are used to 
represent the initial dense region. The chemical composi- 
tions are assumed to be the cosmological residual values 
(IGalli fc Palla Ifl998h . 

We start numerical simulations from this initial setup 
without any radiative feedback effects. Then, the dense 
clump collapses in run-away fashion because of the self- 
gravity and the radiative cooling by H2 . When the cen- 
tral density exceeds a certain limit, n on , we ignite a 12OM 
POPIII star that is located D pc distant from the center of 
the dense clump. The luminosi ty and the effective temper- 
ature of the star are taken from lBaraffe, Heger fc Wooselv I 
(|200lf ). We employ various parameters in the range 
3xl0 2 cm -3 < n on < lxl0 5 cm -3 , and 20pc < D < 200pc. 
In this paper, the ionizing photons are not included, since 
we focus on the photodissociation process. Simulated 
physical time is 10 T yr, which is significantly longer than 
the life time of a 12OM star (~ 2Myr). During the simu- 
lation time, we do not turn-off the source star, as we are in- 
terested in both of the effects by diffuse background/local 
LW radiation (diffuse radiation never turns-off). 

Numerical simulations are mainly carried out on paral- 
lel computer system called "FIRST" installed in Univer- 



3 



sity of Tsukuba. The newly being developed PC cluster 
"FIRST" is designed to elucidate the origin of first gen- 
eration objects in the Universe through large-scale sim- 
ulations. The cluster consists of 240 nodes ( each node 
has dual Xeon processors ). The most striking feature of 
the system is that each node has small GRAPE6 PCI- 
X board called Brade-GRAPE. The Brade-GRAPE has 
four GRAPE6 chips per board (thus 1/8 processors of a 
full board). Th e board has basically the same fun ction 
as GRAPE-6A (jFukushige. Makino fc Kawai Il2005[ ). but 
it has an advantage that it can be installed in 2U server 
unit of PC clusters due to the short height of its heat sink. 

4. NUMERICAL RESULTS 

4.1. Simulation without radiative feedback 

Before we proceed to the results with radiative feedback, 
we give an outline of the results without the feedback pro- 
cess. Since we have no LW flux in these particular runs 
and total mass and the initial density are fixed, the only 
free parameter of these calculations is the initial temper- 
ature, Ti n i. FigJT] shows the snapshots of radial density 
distributions of the collapsing clouds with Ti n i = 100K 
and Ti n i = 350K. The snapshots are taken when the cen- 
tral dense regions are very close to the resolution limits. 
For both of the runs, the distribution roughly converges 
to the p oc r _2 relation, which is we ll known as the Larson- 
Penston type similarity solution ([Larson 1 1 19691 ; iPenstonl 
1969) for isothermal gas. Remark that the larger Ti n i re- 
sults in smaller size of the cloud, since higher temperature 
promotes faster reaction to produce H2, which leads to 
efficient cooling. 

Fig|2] represents the evolution of the cloud core on den- 
sity v.s. temperature plane. Four tracks correspond to 
Ti n i = 50, 100, 200, 350K, respectively. Among these, the 
run with Ti n i = 350K is relatively close to the phys- 
ical states of the collapsing clouds at > 10 3 cm 3 , 
developed from the CDM density fluctuations (see fig- 
ures 11,12 in O'shea & Norman 2006). The runs with 
Tmi = 50, 100K also could be realized in the shock com- 
pressed layers in larger halos with M > 1O 8_9 M0 since 
the gas is cooled down to < 100K due to the efficient 
Ho coolingjShapiro fc Kang II 19871: [Kan^fc Shapiro 1119921 : 
ISusa et al . 1998; Uehara & Inutsuka 20Q0f). Thus, the var- 
ious Tini should be regarded as the variety of the initial 
conditions for the sites of primordial star formation. 

4.2. Typical results with radiative feedback 

First, we show two typical results among various runs. 
Figure [3] shows the density / H2 map on the slice that 
contains the axis of symmetry in the run with T| n i = 
100K, n on = 10 3 cm _3 ,D = 40pc. Top panels represent 
three snapshots, those correspond to t = lOyr, 10 4 yr and 
10 6 yr, respectively. Bottom panels also show yn 2 distri- 
bution at same physical time. In this run, H2 are totally 
photodissociated by the nearby star, even at the center of 
the collapsing cloud (lower panels). H2 fraction never re- 
covers through the simulation time. As a result, the cloud 
cannot keep collapsing because of the absence of coolants. 
The cloud bounces after the adiabatic compression phase. 

On the other hand, Figure [H shows the case same as Fig- 
ure [3] except D = lOOpc. Since the Lyman- Werner band 
flux is smaller than the previous case, H2 is self-shielded 



and survive at the center (lower panels), followed by the 
collapse of the core (upper right panel). Figure [5] shows 
the time evolution of the central H2 fraction for four mod- 
els including above two. The four models correspond to 
D = 40, 60, 80, lOOpc respectively, whereas other parame- 
ters n on = 10 3 cm -3 and Ti n i = 100K are common. The 
horizontal axis shows the time after the collapse started, 
while the vertical axis shows the H2 fraction at the center 
of the cloud. We find that H2 molecules are photodis- 
sociated just after the ignition of the source stars for all 
models. After the dissociation, H2 fraction recovers imme- 
diately. For the models with D = 80, lOOpc, the recoveries 
are sufficient to cool the gas (?/H2 ~ 10 -4 is required), 
whereas it is not true for D = 40pc case. D = 60pc is the 
marginal case, in which the cloud manage to collapse al- 
though it is delayed significantly, because the cooling time 
is slightly longer than the free-fall time. We will discuss 
this delicate behaviour of H2 fraction in section 15.21 

4.3. Summary of numerical runs 

Figure [6] shows the results of various runs. Four panels 
correspond to Ti n i = 50, 100, 200, 350K, respectively. Hori- 
zontal axes represents the turn-on density n on and vertical 
axes show the distance between the source star and the 
cloud center, D. Each symbol in the panels corresponds 
to each run. In the runs denoted by vertices, the clouds 
cannot collapse until the end of the simulation, i.e. 10 7 y r 
after the ignition of the source star. On the other hand, 
the open circles represent the runs with successful collapse. 
Here "collapsed cloud" is defined as the one whose central 
density exceeds the numerical resolution limit determined 
by the Jeans condition. The nominal density correspond 
to the resolution limit is ~ 10 6 — 10 7 cm -3 in the present 
paper. 

The results depend upon the initial temperature Ti n i i.e. 
initial coldness of the collapsing cloud. Initially colder gas 
can collapse for wider range of parameters than those in 
hotter initial conditions. In fact, considering the case with 
Tini = 350K (lower left panel), which is close to the prestel- 
lar core formed directly from CDM density perturbations, 
LW feedback is inefficient for > 10 4 cm -3 for D < 50pc. 
Thus LW feedback from the stars in the same halo would 
not be effective for such dense self-gravitating cores, since 
the virial radius of the ~ 1O 6 M halos at z ~ 20 — 30 is 
several ten parsec. 

On the other hand, considering Ti n i = 50, 100K cases, 
LW feedback is also ineffective (upper panels) even in 
less dense clouds (tin > 10 3 cm -3 ). Thus, if we regard 
these as the stars formed from the shocked layers or fossil 
HII regions, such POP II.5 starsjjh onson fc Bromm 1 12006: 
iGrief fc Bromm 1 12006) would not be affected by LW radi- 
ation even if the source stars are in the same halo. 

These results are discussed physically in detail in the 
next section. 

4.4. Effects of diffuse LW flux on dense cores 

Before we move on to the detailed physical arguments, 
we evaluate the effects of diffuse LW radiation field on the 
primordial star formation. Previous analysis on this issue 
is basically based upon the optically thin approximation 
(e.g. Machacek, Brya n fc Abel 2001), which is only valid 
for low density clouds. lYoshida et al. I (|2Q03h also address 
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this issue including the effects of self-shielding in their cos- 
mological simulations, although in an approximated man- 
ner. Present calculation can give a better understanding 
on this issue especially for dense clumps. 

The intensity of the diffuse LW radiation back- 
ground formed after the first generation stars, is 
expected to be Jlw S, 10~ 21 erg/s/Hz/ c m 2 /str 
for z > 20 (Ciard i. Ferrara fc Abel I l2QQ0t 



iMeisinger, Bryan, fc Haiman II2QQ6D . If we interpret this 
intensity as the flux from distant 12OM single POPIII 
star, the distance D should satisfy D > lkpc. 

In view of the numerical results in Figure El diffuse ra- 
diation field cannot prevent the self-gravitating prestellar 
cores (tin > 10 3 cm -3 ) from collapsing, although the dif- 
fuse field keep irradiating the cloud more than two million 
years. However, it would be possible to stop the collapse at 
much lower density, at which the clouds are not expected 
to be self-gravitating (e.g. Ahn & Shapiro 2006). In such 
phase, dark matter gravitational force dominates the con- 
traction, which is not included in the present paper. It is 
beyond the scope of our present calculations. 

5. ANALYTIC ESTIMATION 

In this section, we describe the analytic collapse criteria 
in the presence of photodissociative feedback by a nearby 
star. The collapse criteria should be t coo \ < i.e. the 
cloud is cooled faster than it shrinks by gravity. After 
some algebra, an analytic expression for t coo \ < t& is de- 
rived, which is compared with numerical results. We also 
compare the derived cooling condition with tdis > te- 

5.1. N on- equilibrium fraction of electrons 

In order to assess the cooling time, we need the num- 
ber density of hydrogen molecules. The amount of hydro- 
gen molecules strongly depends on the electron abundance, 
since they are mainly formed through following reactions: 

H + e-->H _ + ftz/ (1) 

H" + H -> H 2 + e~ 

Thus, we also have to assess the amount of electrons, which 
is out of equilibrium. The fraction of electrons (y e ) at the 
center of the collapsing cloud is not in chemical equilib- 
rium, which follows the following rate equation: 



dye , 2 
— -k rec n N y e . 



(2) 



Here the collisional ionization term is omitted, because the 
temperature of the core is much lower than 10 4 K above 
which the collisional ionization becomes important. Since 
the central part of the core collapses with free-fall time, 
the evolution of the hydrogen nucleon number density 
is described as, 

-1/2 

(3) 



1 dn^ 
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tin dt \32Gmp7iN, 
Combining above two equations, we have 



dy e 
dn^ 



-V2„ ,2 



2/e fc r- 



3tt 



32Gra p ' 



(4) 



k rec is a function of time, because it depends on the tem- 
perature (oc T -1 / 2 ). However, the change of temperature 

2 Remark that equation is valid for N^ 2 > 10 14 cm 2 . 



is not so significant during the collapse, whereas the den- 
sity gets larger by several orders of magnitude. Thus, we 
can approximate the recombination rate k rec as a constant 
when we integrate the equation The equation (jlj) has 
an analytic solution: 

Ve = — " ; " ; / = ( 5 ) 



where y e o denotes the initial electron fraction, and t^no 
denotes the initial number density of hydrogen nucleus. 
Consequently, we obtain the expression of y e as a function 
of number density n^- The asymptotic behaviour of the 
solution ([5]) for ^> ^no is 



1 1 



32Gra D 



k rec V 3tt 

This expression is independent of initial abundance y e o 
and initial density tino, except that y e cannot exceed y e o 
since we neglect the ionization processes. In fact, y e given 
in equation ([5]) are plotted for various initial conditions 
(^no> 2/eo) m Figure H We find rapid convergences of elec- 
tron fractions to equation (|6|) as the collapse of the clouds 
proceed. Therefore, as far as we consider > 10 3 cm -3 , 
we can safely assess the electron fraction by equation (|6]). 

5.2. H2 abundance 

Because of the intense LW radiation, H2 at the center of 
the cloud is basically in chemical equilibrium with given y e . 
Therefore we assess the H2 fraction yn 2 assuming chemical 
equilibrium. Equating the formation rate of H2 and the 
photodissociation rate, we have 

yn 2 = —r (7) 

^2step 

where /c H - denotes the reaction rate of react ion ([1]), festep 
is photodissociation rate by Solomon process. These rates 
are 

(8) 

\ (9) 



k u - = 1.0 x 10" i8 TcmV 

&2step = 1-13 X 10 8 F LW o/sh 



^Vh 2 



lO^m-^ 

where Flwo is the LW flux from the star in the absence 
of shielding, ,f s h denotes the s elf-shielding function derived 
bv iDraine fc Bertoldi I (|l996l ). It is defined as 

, . , _ r 1, x<i 

hh[x) ~ { x~V\ X > 1 

Nn 2 denotes the hydrogen column density of the collaps- 
ing core. Since the core size is approximately given by the 
Jeans length, 7Vh 2 is given as 



^Vh 2 



Aj 

™n2/h 2 — 



1 / 7t/c B T 

2 \/ Gmin^ 



(10) 



Combining equations (|6])- (fT0]h we obtain 

^LWO 



2/H 2 



4.3 x 10" 



2 x 10 _17 erg s _1 cm -2 Hz -1 
7/2 / t \ 11/2 



Vl0 4 cm- 3 y V10 3 K/ 



,10 4 cm- 3 ,' 1 ' if ^ 1 ^ 
Thus we have obtained the formula to assess H2 fraction 
with given LW flux, density and temperatur43 • It is worth 
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noting that y^ 2 is quite sensitive to density and temper- 
ature. If we consider the adiabatic collapse case, the re- 
lation yn 2 oc T 43 / 4 = T 10,75 is satisfied, because of the 
adiabatic relation tin oc T 3 / 2 . We have already found this 
delicate behaviour of yn 2 in the previous section, where 
the very quick recovery of H2 fraction is found in the adi- 
abatic collapse phase (Figure [SJ). Remark that this strong 
dependence on temperature/density holds only if the H2 
is in chemical equilibrium. The relation is not correct for 
the very dense regions in which the LW shielding is signif- 
icant and the dissociation time scale is much longer than 
the free-fall time. However, H2 fraction is large enough to 
cool the gas cloud in such dense regions even without such 
delicate behaviour. 

5.3. Cooling criteria 

Now we are ready to find the cooling criteria. We can 
evaluate the cooling time utilizing the equation ([TT]). The 
cooling function around 10 3 K is approximated as 

( 4.0 x 10- 25 n^H 2 T 3 4 , n N <C 10 4 cm- 3 
Ah 2 ^ < 

[ 2.5 x 10- 21 n N ^H 2 T 3 4 - 7 , n N > 10 4 cm- 3 

(12) 

Here the unit of Ah 2 is erg cm 3 s x , and T3 denotes 
T/10 3 K. This formula is obt ained by fitting the r esults of 
high and low density limit in lGalli fc Palla I (|l998h - Then, 
we can write explicitly the collapse criteria t coo \ < ts as 



< 



3tt 



2Ar 2 (n N , T, F LW ) y 32Gra p n N 



(13) 



Note that Ah 2 is a function of Flw through yn 2 . 

Combining equations ([TT]) . ([T2]) and ([T3]) we obtain the 
cooling criteria as follows: 



T > < 



7.4 x 10 2 



n N 



10 4 cm- 3 



-0.47 



^LWO 



2 x 10- 17 cgs 



0.47 



7.9xl0 2 K ( 4 nN Y 32 ( 
V10 4 cm- 3 / V 



for tin <C 10 4 cm 3 

^LWO 



2 x 10- 17 cgs 



0.43 



for n N > 10 4 cm -3 

(14) 

Therefore, once the temperature at the center of the 
cloud satisfies above condition, the cloud cools and col- 
lapses even in the presence of photodissociative radiation 
field, irrespective of its history on density-temperature 
plane. 

5.4. Comparison of the numerical results to the analytic 

evaluation 

Now the cooling conditions are shown in Figures [8] and 
[9] for two cases on density- temperature plane. Figure [8] 
corresponds to the case in which a 12OM POPIII star is 
located 20 pc distant from the center of the collapsing core, 
whereas it is 80 pc in Figure [9j The horizontal axes denote 
the density at the center of the collapsing core, vertical 
axes are the temperature. The upper right hatched re- 
gions satisfy the cooling condition. The cooling condition 



^cooi < tfi ( corresponds to the hatched region) is obtained 
by directly solving equation ([T3]) . Note that the whole 
cooling region and its boundary satisfy Nn 2 > 10 14 cm 2 , 
which is required to derive equation ([TT]) . Thus, equa- 
tion (p~3|) based on equation ([TT]) correctly describes cool- 
ing condition. 

The dashed curves starting from the left side denote the 
evolutionary tracks on central density-temperature plane, 
with no feedback effects for various initial temperatures. . 
The solid lines which depart from the no feedback curves 
represent the cases where the photodissociation feedback 
is taken into consideration. In fact, the departing points 
from the dashed curves correspond to the timing when the 
nearby POPIII star is ignited. 

We have two distinct classes of runs. First one is the 
case the LW band photons are effectively shielded. Thus, 
the cloud core is hardly affected by the feedback effect. 
Most of the runs in Figure [9] correspond to this case. 

On the other hand, we also have the runs in which the 
LW band photons are not shield enough at the onset of the 
radiative feedback, which is immediately followed by the 
adiabatic collapse (nN oc T 3//2 ). After this adiabatic phase, 
we also have two distinct classes of runs, one of which 
bounces, and the other collapses. The fate of the cloud is 
obviously determined by the cooling condition which we 
derived in this section: if the central density and tem- 
perature of the collapsing cloud can reach the cooling re- 
gion, the collapse continues, otherwise it bounces. In other 
words, in case the cloud can shrink by the inertia until the 
cloud satisfies the cooling condition, the cloud can keep 
collapsing. If not, the collapse is stopped by the enhanced 
thermal pressure during the adiabatic contraction. 

This results basically hold for both of the figures [8] and [9l 
in which different flux is assumed. It is also worth noting 
that the cooling criteria is independent of the initial cloud 
density and temperature, since it is written by some com- 
bination of various reaction rates and physical constants. 
In fact, all the calculations starting from different initial 
conditions evolve following the cooling criteria. 

5.5. Comparison with the condition tdis > tft 

Now we co mpare the present result s with simple criteria 
suggested bv lGlove~ fc Brand I (|200lf ). i.e. t dis > t S - How- 
ever, it is difficult to compare our results directly with 
theirs since they assume fixed uniform cloud density and 
H2 fraction. Therefore, we re-derive the condition for given 
n on and D which is equivalent to tdis < te, utilizing the 
relations found in the previous section. 

Critical condition tdis — is obtained from the assump- 
tion that H2 is in equilibrium: 

&2ste P 2/H 2 = n N y e k u - (15) 
and the dissociation time equals to the free fall time: 



3tt 



^2step - 1/ 32Gm 



(16) 



Multiplying the two equations ([T5]) ([TO]) , and using equa- 
tion (|6|), we obtain 



yn 2 



(17) 



Combining equations ([TT]) , ([9]) and ([TO]) , critical distance 
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that satisfy t^ s = % for a given core density and tem- 
perature T is described as 

DcT = U3pC (l0 24 e7g/s) (lOW 3 ) (gOOk) 

(18) 

Here we use the usual relation between the luminosity 
of the star(I/Lw) and unshielded flux (Flwo)? 

4ttF LW0 D 2 ct = L LW (19) 

Figure [10] shows the loci on n on — D plane along which 
^dis = tfi is satisfied (i.e. D = D cr ) at the ignition time 
of the source star. Two lines correspond to the core 
temperatures (in equation [18]) equal to 300K and 400K. 
If we consider Ti n i = 100K case, the core temperatures 
at the turn-on time of the source star roughly satisfy 
300K <T < 400K (see Figure [2]). Thus, the actual condi- 
tion tdis — tff would be between the two lines. Numerical 
results with T ini = 100K are also superimposed on the 
n on — D plane as the Figure [6] 

Comparing the condition tdis = with the numerical 
results, we find that 1) the simple condition tdi S = is n °t 
so bad, 2) but it overestimates the critical distance by a 
factor of 2-3. For comparison, we also plotted the distances 
that are factor C times smaller than D cr with T = 400K. 
Five lines correspond to C = 0.1,0.2,0.3,0.4,0.5, respec- 
tively. We find C = 0.3 — 0.4 fit the numerical results 
better than the original condition D — D cr does. 

The disagreements between the condition tdis — tfi an d 
the numerical results basically comes from the dynamical 
effects of the collapsing gas. As discussed in the previous 
subsections, H2 fraction in the collapsing core recovers very 
rapidly following the equation (fTT]) during the adiabatic 
compression phase. As a result, dynamically contracting 
clouds can keep collapsing even if condition tdis < tft is 
satisfied when the source star is turned-on. 

Figure [11] also shows the comparison of numerical re- 
sults for Ti n i = 350K with the criteria tdis = te- Since 
the core temperature at the turn-on time is approximately 
100 — 200K for this case, we plot the loci tdis = for 
T = 100K and 200K. In this case, we also find qualita- 
tively the same results as above, however, the disagree- 
ment is slightly smaller than the previous case. In fact, 
we also plotted the loci along which D = CD cr (200K) are 
satisfied for various C. In this case, C = 0.5 — 0.6 gives a 
better fit to the numerical results. The smaller difference 
is due to the fact that the core temperature is so low at 
the turn-on time, that the adiabatic collapse is not enough 
to drive abundant H2 formation. 

In conclusion, tdi s > tff(i.e. D > D cr ) gives a rough 
estimates, but we cannot neglect the dynamical effects on 
the H2 formation. 

6. DISCUSSION 

Present calculations are based upon the assumption that 
the collapsing gas is self-gravitating. This assumption 
would be correct for the prestellar cores formed by some 
fragmentation process of the shocked gas in large halos 
with T v i r > 10 4 K. It is also relevant for the the collaps- 
ing dense prestellar cores (tin > 10 3 cm 3 ) at the center 
of the very first halos formed directly from CDM density 
fluctuations. However, it would not be correct for the ear- 
lier phase (i.e. low density phase) of such prestellar cores 



with nN < 10 3 cm -3 , since dark matter halos are domi- 
nat source of gravity at such early phase. Thus, we have 
to con sider the effect of dark matter for such low density 
cores. lAhn fc Shapiro I (|2006f ) have performed one dimen- 
sional radiation hydrodynamics simulations on this issue. 
They take into account the effects by photoionization as 
well as the LW radiation. They also consider the effects 
of the age of POPIII stars. As a result, they find that net 
effect of radiation by a nearby star is approximately even, 
i.e. negative feedbacks balance with the positive feedbacks. 
These results could also be checked by our code in future. 

LW-band radiation is basically the sum of line emis- 
sions, which should be treated carefully. Unfortunately, 
it is almost impossible to solve line transfer problem di- 
rectly in huge three dimensional simulations by present 
computational resources. Therefore, we have employed the 
self-sh ielding function introduced by iDraine fc Bertoldil 
(1996). However, their formula is basically derived by the 
calculation for static slab, which might be irrelevant to the 
present issue. It is obvious that the self-shielding function 
is inappropriate to represent the absorption of the radia- 
tion in case the motion of the fluid element is supersonic. 
Therefore, the absorption by the gas in the envelope of 
the collapsing cloud would be over estimated , since the in- 
falling envelope is supersonic ([Larson I [l969) . However, if 
the gas motion is not supersonic, the function gives a good 
estimate for the absorption of LW photons. In fact, in our 
numerical simulations, most of the absorption takes place 
in the collapsing core which is subsonic. Thus, we guess 
that we evaluate the absorption of LW radiation properly. 

The effects of ionizing photons are obviously quite im- 
portant, although it is not included in the present pa- 
per. Photoheating effects followed by the photoioniza- 
tion can leads negative feedback effects such as the pho- 
toevaporation of the cloud (e.g. Yoshida et al. 2006), 
creation of shock fr ont (M-type I-front) w hich blow 
out the whole cloud ( Sus a fc Umemural 120061 ), whereas 
it also has p ositive feedback effects such as enhanced 
H2 formation (|Ahn fc Shapiro! l2006f h or formation of 
"H2 barrier" between the source and the collapsing 
cloud (|Susa fc Umemurall2006l ). This is a quite complicated 
problem, however, we also have performed numerical sim- 
ulations on this issue. The details will be discussed in the 
forthcoming paper. 

7. SUMMARY 

We perform radiation hydrodynamics simulations on the 
radiative feedback of POPIII stars. We obtain the well- 
understood collapse criteria of the primordial prestellar 
core both analytically and numerically in the presence of 
Lyman- Werner band photons from nearby stars. The cri- 
teria dictate the importance of self-shielding effects cou- 
pled with hydrodynamics. Consequently, the LW radia- 
tion from a POPIII star cannot halt the collapse of the 
dense clump (tin > 10 3 — 10 4 cm -3 ) even if they are in the 
same halo. We also evaluate the effects of diffuse LW ra- 
diation background, which also is not important for the 
collapse of dense cores with > 10 3 cm -3 . 
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ing and comments. We also thank N. Yoshida and 
M.Umemura for discussions and careful reading of the 
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Fig. 1. — Density distribution of the collapsing cloud without radiative feedback. Two groups of dots represents the runs with Tj n i = 100K 
and 350K, at the epoch when the central part go below the resolution limit. 
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Fig. 2. — Evolution of the central core on density-temperature plane for various initial conditions. In these runs, the radiative feedback 
effects are not included. 
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Fig. 3. — Time evolution of the density distribution and H2 distribution are shown for D = 40pc, n on = 10 3 cm -3 , Ti n i = 100K. Top panels 
represent the density distribution on the slice which include the axis of symmetry after the ignition of nearby star. Three panels correspond 
to t = lOyr, 10 4 yr and 10 6 yr, respectively. Bottom panels also show the snapshots of yn 2 distribution at three epochs. 
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Fig. 4. — Same as Figure [3] except D = lOOpc. 
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Fig. 5. — H2 fraction at the center of the cloud is shown as functions of time. Results from four runs are shown, in which n on = 
10 3 cm -3 , Ti n i = 100K, and D = 40, 60, 80, lOOpc, respectively. The horizontal axis denotes the time after the collapse of the cloud is started. 
The ignition time of the nearby source star is marked by the arrow. Two solid line represent the case in which the cloud fails to collapse 
(D = 40pc) and delayed collapse(Z) = 60pc), whereas the successful collapse cases are shown by the two dotted lines(Z) = 80, lOOpc). 
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Fig. 6. — Summary of numerical runs are shown on n on — D plane. Four panels correspond to the runs with Tj n i = 50K , 100K, 200K and 
350K, respectively. Open circles denote the runs in which the clouds collapse successfully, whereas vertices represents the failed collapse. 
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Fig. 7. — Evolutionary tracks of electron fraction y e are plotted against number density of the collapsing cloud. Each curve corresponds to 
a certain initial condition (n-^Q,y e o). 
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Fig. 8. — Evolutionary tracks on the central density v.s. temperature plane of the collapsing core are shown. Horizontal axis denotes the 
density of the core, whereas the vertical axis denotes the temperature. Dashed lines are the results of runs without radiative feedback, staring 
from four different initial conditions. Solid lines denote the case the results with radiative feedback for various ignition densities n on and 
D = 20pc, which in turn collapse successfully or fail to collapse. In the upper right shaded region, the condition % > t coo \ is satisfied, whereas 
the upper right region bounded by thin dotted line shows the condition where the H2 column density of the core is larger than 10 14 cm — 2 , 
assuming equation JTT). 
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Fig. 9. — Same as Fig. [8] except that D = 80pc. 




Fig. 10. — The conditions t^is = tff(D = D cr ) are plotted by solid lines on n on — D plane for T = 300, 400K. Symbols represent the 
numerical runs. Open circles denote the runs in which the clouds collapse successfully, whereas vertices represents the failed collapse. The 
dotted lines represents the loci along which D = CD cr (400K) are satisfied. C = 0.1,0.2,0.3,0.4,0.5 corresponds to each line. 
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Fig. 11. — Same as Figur JTol except that Tj n i = 350K and the assumed core temperatures are 100K and 200K. The dotted lines represents 
the loci along which D = CD cr (100K) are satisfied. C = 0.3, 0.4, 0.5, 0.6, 0.7 corresponds to each line. 



